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ABSTRACT 

We examine the evolution of the black hole mass - stellar velocity dispersion (Mbh-c) relation over cosmic 
time using simulations of galaxy mergers that include feedback from supermassive black hole growth. Our 
calculations have been shown to reproduce the power-law scaling of the local, redshift zero Mbh-c relation 
and yield the observed normalization for an appropriate choice of the feedback efficiency. In the present paper, 
we consider mergers of galaxies varying the properties of the progenitors to match those expected at various 
cosmic times. In particular, we examine a range in redshifts z = 0-6, where we modify the virial mass, gas 
fraction, interstellar medium equation of state, surface mass density, and concentration of the dark matter 
halos. We find that the slope of the resulting M^-a relation is the same at all redshifts. For the same feedback 
efficiency that reproduces the observed amplitude of the Mbh-c relation at z = 0, there is a weak redshift- 
dependence to the normalization, corresponding to an evolution in the Faber-Jackson relation, which results 
from an increasing velocity dispersion for a given galactic stellar mass. We develop a formalism to connect 
redshift evolution in the Mbh-c relation to the scatter in the local relation at z = 0. For an assumed model for 
the accumulation of black holes with different masses over cosmic time, we show that the scatter in the local 
relation places severe constraints on the redshift evolution of both the normalization and slope of the Mbh-c 
relation. Furthermore, we demonstrate that the cosmic downsizing of the black hole population introduces a 
black hole mass-dependent dispersion in the Mbh-c relation and that the skewness of the distribution about the 
locally observed M^-a relation is sensitive to redshift evolution in the normalization and slope. Improving 
the statistics of the local relation to measure the moments of its distribution will therefore make it possible 
to constrain the evolution of the M^-a relation and black hole population at redshifts that are currently not 
observationally accessible. In principle, these various diagnostics provide a method for differentiating between 
theories for producing the M^-a relation. In agreement with existing constraints, our simulations imply that 
hierarchical structure formation should produce the relation with small intrinsic scatter as the physical origin 
of the Mbh-c enjoys a remarkable resiliency to the redshift-dependent properties of merger progenitors. 

Subject headings: galaxies: formation - galaxies: evolution 



1. INTRODUCTION 

Galactic spheroids exhibit a constitutive relation be- 
tween their stellar kinematics and supermassive black holes 
(SMBHs) through observed correlations betw een black hole 
mass Mbh and either the bulge ma s s Mh„w tMiv oshi et all 
[ 19951 Kormend y & Richstonelll995t lEckart & Genzellll997t 
iFaber et al.ll997HMagorrian et al.ll998UGhez et al.ll998l). or 
the velocity dispersion a {M^—cr, Ferrarese & Merritt 2000; 
Gebh ardt etalll200Ci . The M BH -a relation locally is a tight, 
power-law correlation with logarithmic s lope originally esti- 
mated to be in the range (3 ~ 3.75-4.8 (Ferrarese & Merritt 
1200(1 iGebhardt etalll2000h . Later stu dies improved the ac- 
curacy of the local correlation (e.g. iTremaine et al.l 120021 
hereafter T02) and explored the relation b etween SMBH 
mass and other galaxy properties more fully (IWandell l2002: 
Ferraresdl2 002: Ma rconi & Hundl2003t feernardi et al.ll2003l 
Ba^s^tall2003lK Efforts to extend our knowledge of the con- 
nection between SMBHs and their host galaxies to higher 
redshifts include attempts to measure SMBH properties at 
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z = 4-6 ( Vestergaard 2004) and to quantify the relation be- 
tween SMBHs and galaxy kinematics in distant, active gala x- 
ies (iShields et all2003tlTreu"et all2004tlWalter et al.12 004). 

Observations of active galaxies at z > have yielded 
ambiguous inferences about the nature of the Mbh-c at 
other redshifts, implying both redshift-dependent ( iTreu et alJ 
12004]: [Walter e^dj|2004) and redshift-independent relations 
dShields et al.H2003l) . Knowledge of the Mbh-cj relation over 
cosmological time is crucial to understanding galaxy forma- 
tion if SMBH feedb ack has a significa nt impact on th is pro- 
cess (e.g. Springel et al. 2005a; Robertson et al. 2005). The 
focus of our present work is to characterize the development 
of the Mbh-c relation over cosmic time with a large set 
of hydrodynamic simulations of galaxy mergers that include 
star formation and feedback from the growth of supermassive 
black holes, to use existing data to constrain redshift evolution 
in the Mbh-c relation, and to motivate specific observational 
tests for the M^-a relation at high redshifts using the next 
generation of large telescopes. 

A variety of theoretical models have been proposed to ex- 
plain the Mbh-c relation, based on phys ical mecha nisms such 
as viscous disk accretion ( Burkert & Silk 2001), adiabatic 
black hole growth ( MacMillan & Henrikse r] 120021). gas or 
dark matter collapse lAdams et alJl2001t iBalberg & S hapiro 
120021 Adams et al. 2003) , stellar ca pture by accretion disks 
( Miralda-Escude & Kollmeier 2005), dissipationless merging 
Jciotti & van Albadall2001t iNipoti et al.ll2003l). unregulate d 
gas accretion (Archibald et al. 2002; Kazantzidis et al. 2005), 
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and the self-regulated grow th of black holes by momentum- 
or pressure-driven winds dSilk & Reesl fl9 98i iKinel l2003t 
Murra y et all 120051 iDi Matteo et alJ 12005: Sazon ov et all 
2005JK An important outcome of the these studies has 
been the recognition that the growth of SMBHs may 
have important consequences cosmologically, especially for 
galaxy formation, motivating work linking the hierarchical 
growth of structure with the evolution or demographics of 
SMBHs through studie s of qu a sar evolu tion ( Haiman &^oeb| 
19981 iRichstone et alJ 119981 iFabiar] Il999t iMonaco et alJ 
2000t iKauffmann & Haehnemi2000t lHaehnelt & KauffmanrJ 
2000[ branato et alJ 1200 It iMenou et all 1200 It iMenci et al] 
20031 IDi Matteo et all 120031 12004 IWyithe & Loebl 120031: 
Wvithe1 l2005t iHopkins et al.ll2005alblcldl) and the for mation 



of normal galaxies ( Volonter i et alJl2003t ICox et alJ l2005a: 
Rob ertson et al.l2005l) .~ 

Motivated b y the long-standing view that mergers form 
spheroids (e.g. iToomrell 19771) . IDi Matteo et all (f2005) stud- 
ied the self-regulated growth of black holes in galaxy colli- 
sions and showed that this process can reproduce the M^-a 
relation at z = for progenitors s imilar to local galaxies. 
The work of Di Matteo et al. (2005) served as an important 
early attempt to calculate the effects of black hole growth and 
feedback on the process of galaxy formation by addressing 
the generation of an M^-a relation in self-consistent hy- 
drodynamical simulations. This calculation has served as a 
foundation for proceeding research that has attempted to ex- 
plore the emerging paradigm that supermassive black hole 
formation is a natur al and essential outc ome of hierarchical 
spheroid formation. Spri ngel et al.N2005al) calculated the role 
of AGN feedback in the reddening spheroids formed in merg- 
ers by truncating residual star formation. Robertson et al. 
(2005) explored the effects of AGN feedback in reducing 
the size of spheroids of disk galaxies for ming in the early 
gas-rich environments of high redshifts. ICox et all {2005a) 
calculated the contribution of AGN outflows to the x-ray 
emission from halo gas in ell iptical galaxies. Importantly, 
Hopki ns et al.l (12005 alblddlelil d) calculated a wide range of 
implications of the co-evolution of spheroids and black holes 
for the interpretation of b oth quasar and galaxy observations. 
Hopk ins"et al.l ( 12005 albld d) demonstrated thai modeling the 
obscured growth of black holes during gas-rich galaxy merg- 
ers leads to a luminosity-dependent quasar lifetime that dif- 
fers substantially from standard assumptions for the time- 
dependence of quasar activity. The luminosity-dependent 
quasar lifetime leads naturally to a new interpretation of the 
quasar luminosity function, where low-luminosity quasars 
are primarily massive black holes in sub-Eddington accre- 
tion phases while quasars with luminosities above the break 
in the quasar lumi nosity function a re accreting near the Ed- 
dington rate. Hopki ns et al.l (l2005e) show that the luminosity- 
dependent lifetime and the new interpretation of the quasar 
luminosity function can explain a wide variety of quasar ob- 
servations by deconvolving the growing black hole population 
from the hard X-ray luminosity function, including longer- 
wavelength quasar luminosity functions, the spectrum of the 
X-ray background, and the obscured fraction of quasars as 
a function of quasar luminosity. The model presented in 
Hopki ns et al . (2005e) also leads to a new model for quasar 
clustering ( iLidz eTaH 120051) and predictions for the evolu- 
tion of the faint-end s lope of the quas ar luminosity function 
( IHopkins et alJl2005gl) . Furthermore, Hop kins et all {2005f ) 
show that the co-evolution of spheroids and black holes in 
these hydrodynamical models can account for the evolution 



of the r ed galaxy sequence, bu ildi ng off the previous calc ula- 
tions of Spring el et al.l l l2005al) and lHopkins et alJ ( I2005el) . 

Many of the previously mentioned calculations utilize 
knowledge of the connection between supermassive black 
holes and their host spheroids over cosmological time 
and throughout the process of hierarchical structure forma- 
tion. The calculations presented in the current paper pro- 
vide a detailed hydrodynamical calculation of the nature of 
this spheroid-black hole connection using cosmologically- 
motivated galaxy models at a variety of redshifts, and have 
already served as essential input into many the aforemen- 
tioned calculations. Independently, the detailed calculation 
of the Mbh-c relation with redshift provides an important 
method for testing the premise of spheroidal galaxy formation 
through the hierarchical merging of galaxies against observa- 
tional constraints from the M^w-cr relation observed locally 
and at high-redshifts. These calculations will also serve as a 
basis for future comparisons between theoretical expectations 
for the co-evolution of black holes and spheroi ds with obser- 
vations of ultra-luminous infrar ed galaxies (e.g. Sande rset alJ 
1988; Sanders & Mirabel 1996) and the submillimeter emis- 
sion from galaxies observed with SCUBA (e.g. ISmail et ail 
1997). 

Specifically, in this paper we use the D i Matteo et a Il (l2005l) 
model to explicitly calculate the behavior of the M^-a re- 
lation with cosmic time induced by the redshift-scalings of 
galaxy properties. Using a large set of hydrodynamical simu- 
lations of mergers that include feedback from supernovae and 
accreting SMBHs, the resulting M-rh-u relation is calculated 
for galaxy models that span a large range in virial mass, gas 
fraction, interstellar medium equation of state, surface mass 
density, dark matter concentration, and redshift. 

We determine the allowed redshift evolution in either the 
normalization or slope of the M BH -cr relation consistent with 
constraints from the observed M^-a relation at z = under 
the assumption that black holes of mass Mbh form at a charac- 
teristic epoch Zf- The redshift evolution of the normalization 
in our simulations yields a dispersion in the z = relation con- 
sistent with the observed scatter. We further relate our simu- 
lations to ongoing efforts to measure the M Bli -a relation in 
high-redshift active galactic nuclei (AGN), and discuss future 
observational tests for normal galaxies at high redshifts with 
the next-generation of large aperture telescopes. 

We describe our simulation methodology in |2] report our 
results in ^3] compare with observational constraints in &\ 
discuss our results in |j5] and summarize and conclude in |6] 
For our cosmological parameters, we adopt a flat universe 
with fl m = 0.3, il A = 0.7, il h = 0.04, and h = 0.7. 

2. METHODOLOGY 

We have performed a set of 112 simulations of merg- 
ing galaxies using the e ntropy-conserving formulation 
( Springel & Hernquisl||2002 [) of smoothed p article hydrody- 
namics iLucvl[l977l: iGing old & Monaghan 1977) as imple- 
mented in the GADGET code (Springel et al. 2001; Springel 
2005c). Each prog enitor galaxy i s constr ucted using the 
method described by Springel et al. (2005b), generalized to 
allow for the expected redshift scaling of galaxy properties. 
For a progenitor with virial velocity V™- at redshift z we de- 
termine a virial mass and radius assuming a virial overdensity 
A v i r = 200 times the critical density, using the relations 
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Table 1. 



Models 


Progenitor V vlr [km s ] 


Redshift z 


Gas Fraction / g;ls 


ISM Pressurization ^eos 


Thermal Coupling »?therm 


Local 


80, 115, 160, 226, 320, 500 





0.4, 0.8 


0.25, 1.0 


0.05 


Intermediate--: 


80, 115, 160, 226, 320, 500 


2, 3 


0.4, 0.8 


0.25, 1.0 


0.05 


High-z 


115, 160, 226, 320, 500 


6 


0.4, 0.8 


0.25, 1.0 


0.05 


Low rjtterm 


115, 160, 226, 320, 500 


6 


0.4, 0.8 


0.25, 1.0 


0.025 



Vvi 



For t he Navarro-Frenk- White concentration l lNavarro et alJ 
of eac h progenitor dark matter halo, we follow 



Bullo ck etaTl lEoOl) and adopt 



M„(z = 0) 



-0.13 



(3) 



where MJz = 0) ~ 8 x 1O 12 /i _1 M is the linear collapse 
mass at the present epoch and H(z) is the Hubble parame- 
ter. The co ncentration and virial radius are then used to create 
iHernauisl (119901) profile dark matter hal os using the conver- 
sion described by Snring eTet al.l ll2005bl) . The halos are pop- 
ulated with exponential disks of mass = 0.041M v jj. and 
fractional gas conte nt who se scalelengths are deter- 
mined using the Mo et al. ( 1998) formalism assuming the spe- 
cific angular momentum content /<j equals the disk mass frac- 
tion »Jd and a constant halo spin A = 0.033. The distribution of 
halo spins in dissipationless si mulations is measured to be in- 
dependent of mass and redshift ( Vitvitska et al. 2002) and our 
chosen value for A is near the mode of the distribution. The 
vertical scale heights of the stellar disks are set t o 0.2/?a, simi- 
lar to the Milky Way scaling ( Sieg el et alJ2002l) . The vertical 
scale heights of the gaseous disks are determined via an inte- 
gral constraint on the surface mass density, a self-consistent 
determination of the galaxy potential, and the effective equa- 
tion o f state (EOS) of the m ultiphase interstellar medium 
(ISM, McKee A O s triteai977l ISnrinTd & Hernmiisril2003l 
The model for the ISM utilized in these simulations accounts 
for the multiphase nature of star-forming gas in galaxies by 
volume-averaging the mass-weighted thermodynamic proper- 
ties of cold clouds and a hot, diffuse component on scales not 
numerically resolved. Sr>ring eT& Hernauistl d2003l) present a 
detailed discussion of the implementation and implications of 
this treatment of the ISM, but the primary dynamical effect 
of this multiphase model is to pressurize the ISM as the hot, 
diffuse component increases the effective temperature of the 
star-forming gas. The effective equation of s tate of the gas 
(for a numerical fit, see Robertson et al. 2004) therefore acts 
to stabilize dense gas in galactic disks against|Toomre| (|1964|) 
instabilities caused by self-gravitation ( Sriringel & Hernauist 
l2003URobertson et al.l200 4: Sprin gel et al.l2005bl) . The pres- 
surization of the ISM can be varied using an equ ation of 
state so ftening parameter <7 E os (for details, see Springel et al. 
2005b), which linearly interpolates between isothermal gas 
(<7eos = 0) and a strongly pressurized multiphase ISM (<?eos = 
1). 

Each initial galaxy model consists of 40,000 gas parti- 
cles, 60,000 dark matter particles, and 40,000 stellar par- 
ticles. The gravitational softening of gas and stellar parti- 
cles is set to 100/i _1 (l +z)~ l pc. The progenitors merge on 
prograde-parabolic orbits, with an initial separation R^t = 
5/? v ii/8 and pericentric distance R pa i = 2R,\. Each progeni- 
tor galaxy contains a "sink" particle with a black hole seed 



of mass 10 5 /z 1 Mq allowed to grow through a mode l based 
on Bondi-Hoyle-Lyttl e ton accretion (iHovle & Lvttletonll939t 
Bondi & Hovle 1944; Bondi 1952), implemented as a sub- 
resolution numerical model by Springel et al. (2005b). The 
Bondi accretion rate onto the black hole is determined by the 
surrounding gas properties. To allow for black hole growth 
the BH particle is allowed to probabi listically a ccrete gas 
particles, but as explained in Springel et al. (2005b) this ac- 
creted gas is added to a reservoir from which the Bondi rate 
is smoothly calculated. The BH particles typically accrete 
~ 100 SPH particles during the evolution of the simulations, 
depending on the details of the galaxy models and mergers. 
In this model, a fraction e r = 0.1 of the mass accretion rate 
Mc 2 is radiated and a fraction % erm = 0.05 of this 10% is 
deposited as thermal feedback into the surrounding gas in a 
kernel-weighted manner. The thermal coupling is calibrated 
to reproduce the nor malization of the local M^-a relation 
(Di Matteo et al. 2005), and we note that other prescriptions 
for the energy deposition from blac k hole feedback ca n pro- 
duce similar Mbh-c scalings (e.g.. iMurrav et al]l2005t) . We 
therefore expect our results to be somewhat insensitive to the 
exact choice for how the feedback energy couples to the gas 
as long as the energy is deposited near the black hole. As 
our intent is to calculate the development of the M^w-cr rela- 
tion over time and to avoid assumptions about the presence or 
form of the M^-a relation at higher redshifts, the progeni- 
tor models do not contain bulges but instead yield an M^-a 
relation self-consistently through the merger process. 

We simulate equal-mass mergers of progenitors appropri- 
ate for redshifts z = 0, 2, 3, and 6 with virial velocities in the 
range V wu = 80-500 km s . For each virial velocity, pro- 
genitors with gas fractions / gas = 0.4,0.8 and EOS softenings 
<?eos = 0.25, 1.0 are simulated. The z = 6 runs are also re- 
peated using a thermal coupling constant of 77thei- m = 0.025 to 
gauge the effect of altering the efficiency with which black 
hole feedback heats the surrounding ISM. Table ^ summa- 
rizes our full set of simulations. The current set of simula- 
tions concerns the M^-a relation generated during the pri- 
mary growth phase of the SMBH, which likely requires a ma- 
jor merger to drive the massive inflow of gas to the central- 
most regions of the colliding galaxies. Future calculations 
of possible evolution in the M^w-a relation should ideally 
involve cosmological simulations of galaxy formation that 
naturally account for the mass spectrum of merger progeni- 
tors. Although such simulations are still in development, we 
note that very minor mergers have been show to produce only 
weak star formation in the more massive progenitor for even 
gas-rich systems as the feeble tidal torquing fails to produce 
strong disk in stabilities in the merging galaxies (ICox et al.J 
in preparation). The ability of the most minor mergers to af- 
fect a change in the M^-a relation in the more massive pro- 
genitor galaxies should also be limited. 

After the completion of a merger, the projected half-mass 
effective radius R e and the mass-weighted line-of-sight stellar 
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FIG. 1 . — Relation between black hole mass Mbh and stellar velocity dispersion a produced by the merging of disk galaxy progenitors appropriate for redshift 
z = 0, 2, 3, and 6. These disk mergers produce nearly the same M^u—cr scaling at all redshifts considered, while the simulated relations display a weak redshift 
evolution in their normalization. For normalization evolution of the form a ~ <t,_q the normalization experiences a redshift scaling in the range § = 0.138 

(relative to the simulated z = normalization, dashed lines) to £ = 0.186 (relative to the o bserved z = normalization, dotted lines). At each redshift we vary 
the gas fraction f g and the pressurization of the interstellar medium ^eos ( see Springel et al. 2005a), and we plot results for Mbh and a in progenitor models 
with gf = 0.4, ffFos = 0.25 (blue circles), gf = 0.8,<Jeos = 0.25 (blue diamonds), gf = 0.4, ^eos = 1-0 (red squares), and gf = 0.8,<?eos =1.0 (red triangles). The 
Tremaine et al. 1 2002) best fit to the local Mbh - "" relation (black line) is plotted for reference. Also shown is the Mbh - relation for z = 6 progenitors calculated 
assuming a decreased thermal coupling for black hole feedback (open symbols). Lowering the thermal coupling for black hole feedback by half results in only a 
small increase in the final black hole mass. 



Table 2. Best Fit M Bii -<r Relations 



Progenitor Redshift 


a 


P 


Alogyv/ BH 


Q /3=4.02 


A /3=4.02 
log jW B h 


z = 


8.01 


3.87 


0.26 


8.01 


0.26 


z = 2 


7.83 


4.10 


0.28 


7.83 


0.28 


z = 3 


7.72 


4.02 


0.26 


7.72 


0.26 


z = 6 


7.44 


3.62 


0.24 


7.45 


0.26 



velocity dispersion a measured within an aperture of radius 
R e are calculated for each remnant galaxy for many viewing 
angles. This determination of a for simulated systems is com- 
mensu rate with the observational method of Gebhar dt et alJ 
(2000) that measures a luminosity-weighted velocity disper- 
sion within an effective radius. The black hole mass M BH 
in each galaxy is tracked during the simulation and com- 
pared with the average dispersion a to examine any resulting 
Mbh— (j correlation. 

3. RESULTS 



Figure^shows the Mbh-t relation produced by the merg- 
ing of progenitor galaxies appropriate for redshifts z = 0, 2, 
3, and 6 (upper left to bottom right). At each redshift, we 
perform a least-squares fit to the Mbh-c relation of the form 

logM BH = a + ^log(cr/cr ) (4) 

where the relation is defined relative to <7o = 200 km s" 1 . The 
observed relation has a normalization coefficient a = 8 . 1 3 and 
slope f3 = 4.02 (T02, solid line). Table lists the best fit a 
and (3 from the simulations at each redshift, along with the 
dispersion Ai og M BH about each best fit relation. Table |2] also 
reports the best fit normalization or =4 02 and resultant disper- 
sion A(^^ assuming the slope of the Mbh-c relation is a 
constant j3 = 4.02 with redshift. The intrinsic Mbh-c disper- 
sion at any given redshift induced by varying the gas fraction 
or ISM equation of state is similar to the observed dispersion 
Aio g M BH = 0.25-0.3 (T02). As Figure Hand Table demon- 
strate, an Mbh-c relation consistent with a power-law scaling 
of index (3 ~ 4 is predicted by our modeling for normal galax- 
ies out to high redshift. We infer that the physics that sets the 
scaling of the Mbh~ct relation is insensitive to an extremely 
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wide range of galaxy properties. Virial velocities of the pro- 
genitors range from V vir = 80 km s" 1 to Vyk = 500 km s" 1 at 
each redshift while the dark matter concentrations range from 
Cvh — 15.6 in the smallest progenitor at z = to C wu ~ 1.5 in 
the largest progenitor at z = 6. For the M^-a relation at each 
redshift, models with moderate ISM pressurization typically 
experience slightly more black hole growth than models with 
strong ISM pressurization. The increased thermal feedback 
in the model with a stiffer equation of state improves the dy- 
namical stability of the gas, lessens gas angu lar momentum 
loss during the merger ( Robertson et al. 200 4 120051). and r e- 
sults in less fueling of the black hole (Hopkins et al. 2005e). 

3.1. The Normalization of the Mbh - c Relation with Redshift 

The normalization of the M^-a relation evolves weakly 
with redshift relative to the observed z = result. To charac- 
terize this, we generalize the parameterization of the M^-a 
relation to the form 

logM BH = a + /31og(cr/cr )-£log(l+z) (5) 

with £ defined to be positive if the characteristic black hole 
mass decreases with redshift or the characteristic velocity dis- 
persion increases with redshift. We determine the value of £ 
for our simulations assuming either the observed z = nor- 
malization a bs = 8.13 (Figure^ dotted line) or the calculated 
z = normalization a ca i c = 8.01 (Figure [2 dashed line) with 
the slope fixed at (3 = 4.02. We find evidence for a weak evo- 
lution, with £ = 0.186 for a \, s and £ = 0.138 for a ca i c using 
the best-fit parameters reported in Table |2] We now consider 
whether the weak evolution we find in the normalization of 
the Mbh-c relation represents a physical effect caused by the 
redshift-dependent properties of the progenitors or is an arti- 
fact of our modeling. 

3.1.1. Systematic Considerations 

The systematic effects of our modeling include the redshift- 
scaling of halo virial properties and an assumed fixed ther- 
mal coupling of black hole feedback to the gas in our simula- 
tions. The form of the progenitor galaxies is inferred from the 
linear theory of dissipational galaxy formation and is there- 
fore idealized, but in the absence of cosmological simula- 
tions with comparable resolution that include feedback from 
SMBH growth, our choice appears well-motivated. While 
the model considered here assumes a fixed thermal coupling 
for black hole feedback, if the thermal coupling changes with 
some bulk property of galaxies with redshift (metallicity, for 
example) then this assumption may influence our results. The 
magnitude of thermal coupling ^ t herm = 0.05 of black hole 
feedback to the surrounding gas is calibrated to reproduce 
the locally observed relation. To illustrate the impact of this 
choice, Figure shows the change in normalization of the 
Mbh-c relation for z = 6 progenitors introduced by lowering 
the thermal coupling parameter to ?7therm = 0.025 (open sym- 
bols, lower right panel). As the coupling of the feedback to 
the gas of the galaxy is decreased, the rate of energy input to 
the surrounding material by a black hole of a given mass and 
accretion rate is lessened. The transition between Eddington- 
limited and self-regulated growth then occurs at a correspond- 
ingly higher black hole mass and, as a result, the final mass of 
the SMBH increases. Decreasing the thermal coupling also 
mildly flattens the slope of the relation for the highest mass 
black holes. Overall, the change in the thermal coupling does 
not re-normalize the M^-a relation for z = 6 progenitors to 




FIG. 2. — Relation between total stellar mass M„ and black hole mass A/rh 
for merger progenitors appropriate for redshifts z = (black circles), z = 2 
(red diamonds), z = 3 (blue triangles), and z = 6 (green squares). The ratio of 
black hole to stellar mass remains roughly constant for progenitors at all red- 
shifts we consider, indicating that the normalization evolution of the Mrh - o" 
relation measured in the simulations is not due to a decoupling of bulge and 
black hole growth. At each circular velocity and redshift, four models for 
the gas fraction and interstellar medium equation of state are considered. We 
plot bulge-dominated remnants only. The MRH - M,ulee relations observed by 
Marconi & Hunl i^f)03) (solid line) and Maeorrian et al.l l!998t) (dotted line) 
are shown for comparison. 



coincide with the T02 result. We also note that varying the 
gas fraction in the progenitors from / gas = 0.4 to / gas = 0.8 has 
little effect on either the normalization or the scaling of the 
relation at any given redshift. 

A sensible measure to gauge a possible evolution of the 
black hole growth with redshift is the comparison of black 
hole mass Mbh with total stellar mass M*, which serves as a 
convenient proxy for the observed MBH-A^buige relation. Fig- 
ure 13 shows the Mbh - M+ relatio n for mergers at a ll red- 
shifts w e consider, with the Marconi & Hunt! ||2P03) (solid 
line) and Magor rian et alJdl99 8) (dotted line) M B11 -M hu \ g(: re- 
lations plotted for comparison. We show the largest stellar 
mass systems, which a re mostly spheroidal and lie near the 
iMarconi & Huntl l )200 3) relation. Any redshift-dependent im- 
pact on the black hole mass at a given stellar mass is not dis- 
cernibly stronger than the black hole mass increase for de- 
creasing ISM pressurization and increasing gas fraction. We 
therefore conclude that the systematic uncertainties of our 
modeling do not generate an unphysical evolution in the black 
hole mass. Instead, it appears more likely that the variation of 
the normalization of the Mbh-c relation with redshift has a 
physical origin in the structural changes of the merging galax- 
ies, a possibility we examine next. 

3.1.2. Physical Considerations 

The redshift-dependent structure of progenitor systems ex- 
pected from the cosmological scalings of galaxy properties 
provides a plausible explanation for the weak evolution in the 
normalization of the Mbh-c relation with redshift. For ex- 
ample, the depth of the potential well in the center of each 
remnant will influence the velocity dispersion, while for the 
regions where the baryons dominate the stellar mass of the 
spheroid determines the potential. The characteristic densi- 
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FIG. 3. — Relation between stellar velocity dispersion a and total stel- 
lar mass M* of remnants for merger progenitors appropriate for redshifts 
z = (black circles), z = 2 (red diamonds), z = 3 (blue triangles), and z = 6 
(green squares). The simulated remnants display an increasing velocity dis- 
persion for a given s tellar mass with redshif t, which may be viewed as a 
weak evolution in the Faber & Jackson 1 1976) relation. The evolution in ve- 
locity dispersion for a given stellar mass induces a corresponding evolution 
in the calculated Mbh _ c relation (see Figure Q. At each circular veloc- 
ity and redshift, four models for the gas fraction and interstellar medium 
equation of state are considered. A total ste llar-mass version of th e Faber- 
Jac kson relation constructe d by combining the Tremaine et al. 1 2002 ) Mbh - a 
and Marconi & Hunt 1 2003) MBH _ ^bulge relations is shown for comparison 
(solid line). 



ties of the progenitors and remnants increase with redshift and 
any relation involving the stellar mass may evolve cosmolog- 
ically. Similar trends can be measured with our simulations. 
Figure [3] shows the stellar velocity dispersion as a function 
of the total stellar mass of each remnant for the simulations 
at each redshift, along with a stellar mass formulation of the 
Faber-Jackson luminosity-stellar ve locity dispersion r elation 
found from the combination of the M arconi & Huntl ( 12001 
MeH-A^buige and T02 Mbh-c relations (solid line). For a 
given stellar mass, the velocity dispersion increases for higher 
redshift progenitors, with the greatest increase for the most 
spheroid-dominated systems. We therefore ascribe the weak 
evolution in the normalization of the Mbh-c relation with 
redshift to a systematically increasing a for a given Mbh or 
M*, and argue that the steeper potential wells of high redshift 
galaxies is a primary cause for this apparent redshift evolution 
in the stellar mass formulation of the Faber-Jackson relation. 

The evolution of the stellar velocity dispersion for a given 
stellar mass between redshifts z = 2 and z = 6 can be com- 
pared with the redshift-dependent properties of our progeni- 
tor models. Changes in the structure of galaxies can lead to 
an evolution in the stellar velocity dispersion that scales with 
redshift in a similar fashion to the evolving normalization in 
our calculated Mbh-c relations at redshifts z = 2-6. A pos- 
sible origin for these effects can be illustrated by the use of a 
simple galaxy model that demonstrates the impact of redshift- 
dependent structural properties on the central stellar kinemat- 
ics. To demonstrate this, we consider the velocity dispersion 
of a stellar spheroid residing within a dark matter halo with 
fixed virial mass M W1I as a function of redshift. For simplicity, 
we assume that the stellar mass within the halo is a fixed frac- 



tion of the virial mass and that the central potential is domi- 
nated by the baryons. We model both the dark matter halo and 
stellar spheroid as Hernquist profiles with the stellar spheroid 
scalelength a s a fraction b of the halo scalelength a^. Equa- 
tions ([0|3j describe the redshift scalings of the dark matter 
halo s for a given virial velo city VW. Incorporating the results 
from lSpringel et afl d2005bl) . we find the spheroid scalelength 
will depend on virial mass and redshift as 

fl s (M vir ,z) = — J 2 [ln(l+C vn -)-C vir /(l+C vir )] (6) 

Cyir 

x[GM vir ] 1/3 [10H(z)T 2/3 

where C v ; r is the mass- and redshift-dependent concentration 
given by Equation Q. The stellar velocity dispersion of the 



Hernquist profile depends on the scalelength as a oc a s 



1/2 



The power-law scaling of a oc (1 +zp for 2 < z < 6 then 
ranges from £ ps 0. 17 for aM vir - 1O 1O M halo to £ « 0.28 for 
M v j r ~ 10 14 Mq, which encompasses the approximate scaling 
£ found in the simulations. Below redshift z = 2, the redshift 
evolution of the stellar velocity dispersion begins to flatten 
and is roughly constant between redshifts z = and z = 1 . The 
relative constancy of the stellar velocity dispersion in this red- 
shift regime is caused by the growing influence of the cosmo- 
logical constant on the Hubble parameter as it balances the 
increase in the virial concentration. We therefore expect the 
normalization of the Mbh-c relation to remain roughly fixed 
for purely dissipationless systems between z = and z = 1 , and 
begin to evolve roughly as a weak power-law in 1 + z toward 
higher redshifts. For the dissipational systems we simulate, an 
additional increase in the stellar velocity dispersion between 
redshifts z = and z = 2 beyond the estimate for purely stel- 
lar systems is needed to fully explain the evolution measured 
in the simulations by redshift z = 2. We leave such an explo- 
ration for future work, but note here that either gas dissipation 
or feedback effects from the supermassive black hole may in- 
fluence the scaling of the stellar dispersion with redshift or its 
mass dependence. 

4. COMPARISON WITH OBSERVATIONS 

4. 1 . Scatter Induced in the Local Mbh _ c 
Relation by Redshift Evolution 

If the Mbh-c relation evolves with redshift, any redshift 
evolution in the cosmic black hole population will affect 
the statistical distribution of galaxies about the observed 
M rh-cx relat i on at z = 0. Following the formalism presented 
bv I Yu & Lul J2004 YL04), we can statistically characterize 
the connection between active galactic nuclei and the local 
Mmi-a relation by assuming BH growth occurs primarily 
during periods of quasar activity. If we further assume that 
the Mbh-c relation generated by mergers at each redshift has 
an intrinsically small scatter, we can calculate the allowed red- 
shift evolution in a given Mbh consistent with the dispersion 
measured in the local Mbh-c relation. 

From logMBH), the probability distribution function 
(PDF) for a black hole of mass Mbh to have formed at a red- 
shift Zf, we can calculate the expectation value of a function 
q(zf,logM BH ) given log Mbh as 

(^|logM B H> = J q(z f ,logM Bll )P(zf\\ogM BH )dzf (7) 

and the moments about (q\ log Mbh) as 

V? n = I [ ? (z / ,logM B H)-(?|logM B H}]" (8) 
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FIG. 4. — Integrated probability distribution function P(z/ > z| logMBH) 
for a supermassive black hole of mass Mrh t o have formed at a redshift 
Zf > z, as determined by the best fit model of Hopki ns et all I2005el) (see 
ifrl.ll for details). Shown is the P{z/ > z\ logMen) for logMRH = 6.0 (solid 
black), 7.0 (purple dotted), 8.0 (blue dashed), 9.0 (green dash-dotted), and 
10.0 (red dash-double dotted). More massive black holes are typically 
formed at successively higher redshifts, reflecting the cosmic downsizing 
trend inferred from evolution in the quasar luminosity function. This P(z/ > 
z| logMRH) distribution is combined with functional forms for redshift- 
dependent Mrh - "" relations to estimate constraints on possible evolution in 
the Mrh - "" relation by comparing with observations at z = (see Figures l5l 
and|SJ. 



xP(z/|logM BH )dz/. 

We can identify the expectation value (cr| logMBH} with the 
Msn-a relation observed at z = 0, and associate the sec- 
ond moment with the measured M BH -cr relation dispersion 
A °o b g S M BH = \/74 = 0.25-0.3 dex (T02). If the velocity dis- 
persion depends on some larger set of independent variables 
x as cr(logMBH,x), and P(zf | log Mbh) is independent of x, we 
can extend this definition of the expectation value to calculate, 
e.g., (<j| logMBH,x). The analysis of YL04 treats a more gen- 
eral scenario where joint PDFs for properties associated with 
black hole mass are functions of more than one variable, and 
we restrict our attention to a single PDF directly relating Zf 
and logMBH- Equivalently, that P(zj | logMBH) is independent 
of x reflects our operating assumptions including: 1) all scat- 
ter in the M B h-<t at z = owes to evolution in a with redshift 
(i.e. if x is null, u does not depend on z and the dispersion 
in the Mbh-c relation would be Ai og M B „ = 0) and 2) galac- 
tic evolution subsequent to the primary growth phase of the 
black hole does not affect the direct connection between z/ 
and log Mbh- 

The Mbh-c relation generated by the merging of progeni- 
tor galaxies appropriate for a range of redshifts demonstrates 
an evolution in its normalization with respect to the observed 
relation at z= 0. By utilizing a known P{zf\\ogMsu) and a 
functional form for er(z, logMBH, x) that captures this evolu- 
tion, the dispersion about the average relation (a\ logMBH, x) 
can be calculated and compared with the observed A°^ Mmi . 
The small observed scatter in the Mbh-c relation will then 
constrain x for a known P(z/\ logMBH)- To this end, we con- 
sider Mbh-c relations with an evolution in normalization of 
the form of Equation (|5j with the identification of x = £, or an 



evolution in the slope of the form 

logM BH = a + /3(l+z) log(cr/a o ) (9) 

with x = <fi. We proceed to calculate the resulting variance {1% 
about the mean Mbh-c relation at z = 0, using these forms for 
o-(M BH ,z) and the observed values of a = 4.0, (3 = 8.13 (T02). 
As the sign of either <f> or £ will strongly affect the symmetry 
of the inferred distribution about the mean relation, we also 
characterize the skewness 71 = /if /(/if) 3 / 2 . 

We now adopt a form for P(z/\ logMBH), using the results 
of lHop kins et al. (2005e, H05), and calculate the probability 
P(z/\ logM B n) for a black hole of mass Mbh to have formed 
at redshift zj- The H05 modeling uses luminosity-dependent 
quasar lifetimes and obscuring column density distributions 
measured from the simulations presented in this paper to in- 
fer the black hole mass function «(Mbh,z) from a host of 
constraints from observations of quasars. The evolution of 
«(Mbh , z) with redshift is primarily determined by an evolving 
characteristic luminosity L*(z) that determines the location of 
the maximum of n{M^(L),z)- In principle, logMBH) 
is well determined in the context of the H05 model for the 
cosmic black hole population via the observed evolution of 
quasar luminosity functions in multiple wavebands and con- 
straints on the mass density of black holes at z = 0. In prac- 
tice, the somewhat weak observational constraint on the evo- 
lution of the break in the luminosity function at the high- 
est observationally-accessible redshifts requires a choice for 
the functional form of the evolving characteristic luminosity 
L+(z) that determines the peak of the black hole mass func- 
tion zi(Mbh(£),z)- As an Ansatz we adopt the form for L*(z) 
chosen by H05 to best fit the available constraints from ob- 
servations, though we note that the calculations presented in 
section §6 of H05 and below can be applied to more general 
forms of «(Mbh,z) and L*(z) and the analysis presented here 
could be repeated for other interpretations of the quasar lumi- 
nosity function than that adopted by H05. 

The resulting P(z/\ log Mbh), plotted in an integral form in 
Figure0]and reproduced from H05, reflects the best fit model 
for the evolution of the cosmic black hole population given the 
available data and the H05 interpretation of the quasar lumi- 
nosity function. Cosmic downsizing, the concept that black 
holes of la rger mass formed c haracteristically at higher red- 
shifts (e.g. Cowie et al. 2003; Steffenetal. 2003), is appar- 
ent in our P(Zf\ logM B n) as the formation-redshift PDF for 
low mass black holes increases strongly at low redshifts. The 
characteristic downsizing of black hole masses as the universe 
ages leads to an immediate general conclusion about possible 
evolution in the M^-a relation: if black holes of differing 
masses form at different characteristic redshifts, any evolu- 
tion in the Mbh-c relation will introduce an MfiH-dependent 
dispersion in the observed relation at z = 0. While subse- 
quent processes may reduce or amplify the mass dependence 
of Ai og M BH , if observations reveal a trend in Ai oe M BH between 
differing mass bins the result could be associated with an evo- 
lution in the Mbh-c relation. Below, we characterize the mass 
dependence of the statistical moments of the distribution of 
galaxies around the mean M^-a relation and relate them to 
possible forms for M^-a evolution. 

4.1.1. Normalization Evolution 

For a normalization evolution of the form of Equation (0, 
the dispersion Ai og M BH an d skewness 71 for the M^-a rela- 
tion at z = can be readily calculated using Equations CHE) 
given a P(z/\ log Mbh) describing the evolution of the cosmic 
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FIG. 5. — Dispersion Ai og M BH and skewness 71 about the local Mbh - & relation caused by evolution in the normalization of the Mbh _ a relation with redshift of 
the form log(MBH) = £*+/31og(o"/°"o)— £ log(l +z). Redshift-evolution in the normalization induces a dispersion in the M^u—tr relation at z = since black holes of 
different masses form at characteristically different cosmological times (see Figure|4}. The dispersion induced by normalization evolution with redshift-scalings 
of |£| = 0.05 (solid), 0.10 (dotted), 0.20 (dashed), and 0.25 (dash-dotted) is plotted. The observed dispersion Aft^ = 0.25 — 0.3 (shaded region) provides a 
tight constraint on the evolution in the normalization, which is limited to approximately § < 0.2. If the normalization evolves, the sign of the skewness in the 
distribution about the local Mbh - °" relation indicates the sign of £. 



black hole population. Figure|5]shows the dispersion in black 
hole mass A] og M BH as a function of SMBH mass introduced in 
the M Bli -a relation at z = by normalization evolution over 
a range in the power law index f (left panel). As Ai ogMBH is 
an even function of f , we plot Ai og M BH f° r different values 
of |f |. The observed range of Af b g s MBH = 0.25-0.3 dex (red 
shaded region) can be compared to the normalization evolu- 
tion measured in the simulations of f « 0.138-0.186 (blue 
shaded region), and the simulations agree well with the ob- 
servational constraints. Given the observational data and our 
chosenf(zy I logMBH), evolution in the normalization with ap- 
proximately |f I < 0.2 is allowed. The mass-dependence of 
AiogM BH directly reflects changes in the black hole popula- 
tion during the evolution of the M^-a normalization with 
redshift, with more massive black holes being formed during 
characteristically earlier epochs when the M-rh-u normaliza- 
tion was more discrepant from the z = relation than dur- 
ing the formation epoch of smaller mass black holes. As 
the most massive black holes form over a wider range in 
redshift than smaller black holes, their broader PDFs also 
contribute to the scatter about the mean relation. If fur- 
ther sources of scatter in the Mbh-c relation are unimpor- 
tant and subsequent galactic evolution does not strongly ef- 
fect P(zj -| logMBH), then the slope of Ai og M BH provides direct 
information about the typical redshift of formation for black 
holes of a given mass. We note that since the black hole and 
spheroids form contemporaneously in our model, the evolu- 
tion owing to f in the calculations also suggests a correlation 
between the deviation of a galaxy from the Mbh-c relation 
and the age of the spheroidal stellar population. For normal- 
ization evolution such a correlation results from the mono- 
tonic mapping of the measured velocity dispersion deviation 
S loga = [logo-(M BH ,z/)-logcr(MBH,z = 0)] = flog(l+z/) to 



a formation redshift Zf, which means the probability distri- 
bution function P((5i oe(T |logMBH) maps monotonically onto 
P(z f \logM BH ). 

Further information on f can be gleaned from the skewness 
71 of the distribution (Figure [5] right panel). For normaliza- 
tion evolution of the form of Equation (0, the resulting skew- 
ness has the same sign as the power law index and immedi- 
ately indicates whether the M B h-c normalization increases or 
decreases with redshift. All positive values of f have the same 
skewness, with the same functional shape as the skewness de- 
fined by negative values of f reflected about the origin. 

4.1.2. Slope Evolution 

While our simulated M^-a relations do not show any 
monotonic trend in their slope with redshift, many theories 
for the Mbh-c relation have some functional depend ence in 
the slope with redsh ift or black hole mass (e.g. S azonov et alJ 
200.4 lWvirhdl2005l) . The disper sion A] gjv/ BH and skewness 
71 for the Mbh-c relation owing to slope evolution of the 
form of Equation l|9} is also accessible using the method 
described for characterizing normalization evolution. Fig- 
ure|6]shows the dispersion Ai og M BH (left panel) about the lo- 
cal Mbh-c relation induced by slope evolution, plotted for 
4> = 0.2 (solid), 0.3 (dashed), -0.2 (dotted), and -0.3 (dash- 
dotted). To illustrate the mass dependence of Ai og M BH an d 7i 
above and below the pivot we choose the fixed intercept to be 
logMBH.fixed = a = 8.13, but the same analysis could be per- 
formed for any intercept. For consistency with the observed 
dispersion A°^ g s MBH = 0.25-0.3 (shaded region), the evolution 
in the slope is constrained to approximately \<f>\ < 0.3 for our 
chosen intercept. If the fixed pivot intercept for a slope evo- 
lution of the Mbh-c relation is small (logMsH ~ 6), then the 
constraint on \(f>\ would be much tighter as Ai og M BH is a strong 
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FIG. 6. — Dispersion Ai og M BH (left panel) and skewness 71 (right panel) about the local Mbh - "" relation induced by an evolution in the slope of the Mbh - "" 
relation with redshift of the form log(MBH) = et+/3(l +z)^ log(<r/cTo). Redshift-evolution in the slope induces a dispersion in the Mbh _ c relation at z = since 
black holes of different masses form at characteristically different cosmological times (see Figure|4J. The dispersion is plotted for = 0.2 (solid), 0.3 (dashed), 
-0.2 (dotted), and -0.3 (dash-dotted). To illustrate the mass dependence of Ai og M BH and 71 above and below the pivot, we choose the fixed intercept to be 
logMflH, fixed = a = 8.13. The observed dispersion Ajj^ = 0.25—0.3 (shaded region) constrains evolution in the slope to approximately \<j>\ < 0.3. If the slope 
of the Mbu-<t relation evolves and the intercept falls in the observable range, the sign of <f> could be constrained by examining the sign of the observed skewness 
as a function of black hole mass because 71 changes sign at the pivot mass. 



function of Mbh above the pivot mass. The skewness induced 
by a slope evolution (Figure |6j right panel) contains more 
information than for normalization evolution. The skewness 
provides a constraint on the pivot mass if the pivot falls in the 
observationally accessible range, as the sign of the skewness 
changes at this intercept in the redshift-dependent relation. 

In principle, any model that produces a redshift evolution 
in the Mbh~ct relation can be compared with the observa- 
tions in a similar manner, and given a P(z/\ logMeH) for the 
cosmic black hole population, distinct models for generating 
the Mbh-c relation could be observationally differentiated. 
We note here that comparisons with observations are most 
meaningful if the black holes in the measured Mbh-c rela- 
tion represent a fair sample of the total black hole population 
for a given mass bin. For instance, an increase in the number 
of bulges in spiral systems with well-determined Mbh and a 
would be helpful. Unfortunately, dust effects make stellar dy- 
namical measurements late-type galaxies difficult and a more 
effective approach for determining black hol e masses in thes e 
systems will need further development (e.g. iHo et alJ l2Q02). 
We encourage further observational efforts to eliminate selec- 
tion biases in the sample used to determine the local Mbh-o" 
relation. 

4.2. Observations ofM^w -1 ^ in AGN at z > 

Our simulations can be directly related to the Mbh-ct rela- 
tion in normal galaxies at redshifts z > 0. However, previous 
observational efforts to measure the Mbh-c relation in galax- 
ies beyond the locally accessible sample have concentrated 
mostly on using the Mrh-ct relation indicators in active galac- 
tic nuclei (Shields et al. 2003: iTreu et al.ll2004j). com bining 
H/3 li newidth estimates for black hole mass (e.g. Wand el et al.l 
1999) and [Olll] linewidths as a proxy for stellar velocity dis- 



persion (e.g. lNelsorl2000a|) or dispersion measurements from 
stellar spectral features (ITreu et all2 001). Other estimates for 
the Mrh-<7 relatio n in high-z quasars have been made (e.g. 
Walter etal]l2004l) by comparing a black hole mass inferred 
from Eddington-limited accretion onto a bright quasar with 
the dynamical mass estimated from Very Large Array obser- 
v ations of mole c ular g as kinematics. 

IShields et all (120031) found that AGN between z = 0- 3.5 
obey an H/3-[Olll] relation consistent with the local Mbh-c 
relation determined by T02, and the observations at any sin- 
gle redshift were also consistent with the z = result, though 
the statistics are somewhat limited. The observations of 
ITreu et all J2004I) imply a conflicting result with seven AGN 
at z ~ 0.37 inferred to primarily lie above the z = Mbh-ct re- 
lation, t hough with s ignificant scatter and large estimated er- 
ror. The Walter et al. (2004) CO(3-2) observations also place 
the z = 6.42 quasar above the z = Mbh-c relation and claim 
to be inconsistent with the presence of a stellar bulge of the 
size needed to match the local MBH-Mb u i ae relations. 

The increased scatter observed in the AGN samples and the 
deviation of high-z quasars relative to our calculations may 
originate in the dynamica l state of the ob served systems (for a 
related discussion, see [yestergaard 2004). If these systems 
correspond to merging galaxies, they may not be dynami- 
cally relaxed and these samples may therefore display a larger 
range of velocity dispersions than do the dynamically relaxed 
systems that primarily comprise the locally observed Mbh-c 
relation or the redshift-dependent Mbh-c relation in relaxed 
galaxies produced by our simulations. For instance, obser- 
vations of local Seyfert galaxies show that tidally-distorted 
hosts have systematically high er emission line widths at a 
given stellar velocity dispersion ( N elson & Whittle! 19961) . An 
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exploration of the M^-a relation during the peak AGN 
phases of our simulations is planned in future work, but given 
the current observational and theoretical uncertainties we do 
not believe our simulations are in conflict with the observa- 
tions. However, our simulations are more consistent with the 
Shiel ds"e{ alj (120 03) claim of no evolution than with either 
the iTreu et al.l (120041) or IWalter et all ( 120041) claims of evi- 
dence for evolution in the M^-a relation. In our model 
where the M^-a scaling is set during spheroid formation, 
the evolution of systems tends to occur either quickly in the 
Mbh direction during the exponential growth phase of the 
black hole or more gradually along the M^w-cr relation dur- 
ing self-regulated black hole growth. The presence of massive 
black holes at high redshift without massive bulges would 
be rare in the context of our modeling, and such outliers 
would likely have to undergo purely stellar major mergers 
to reach the loca l Mbh-ct relation. We also note that if the 
ITreu et all J2004I) offset of S log a = -0.16 measured for AGN 
at z = 0.37 holds for the M^-a relation generated by the 
merging of normal galaxies at that redshift, then a power-law 
evolution in the normalization of a would imply £ = -1.17. If 
some subsequent process does not reduce the deviation from 
the z = Mbh-c relation for a given galaxy, then with our 
P(zfllogMpH) an d the evolution £ = -1.17 inferred from the 
ITreu et all (120041) measurements, the scatter in the Mbh-c re- 
lation at z = would be Ai ogMBH = 0.92-2.215 over the range 
logM BH = 6-10. This scatter is much larger than both the 
observed dispersion and the dispersion predicted by our sim- 
ulations and formation redshift distribution P(zj | logM B H). 

4.3. Observations of Inactive Galaxies at z > 

For the next-generation of large telescopes, a measurement 
of the central stellar velocity dispersion in normal galaxies at 
high redshifts may be feasible, and our simulations can con- 
nect directly with such observations. A proper census of the 
cosmic black hole population as a function of redshift will 
necessarily involve an estimate of the mass of SMBHs at the 
centers of inactive galaxies. Without black hole activity to 
provide an estimate of the black hole mass, and given the 
likely prohibitively extreme resolution (< 10 pc) needed for 
stellar dynamical estimates of SMBH masses, inferring black 
hole masses through the Mbh-c relation via stellar velocity 
dispersion measures for the cosmological spheroid population 
may provide the best way of constraining the evolution of the 
black hole population with redshift. These measurements will 
necessarily involve an assumed model to relate Mbh and a as 
a function of redshift, such as the relation calculated by our 
simulations. Our work then provides additional theoretical 
motivation for attempting these observations with future gen- 
erations of large infrared-capable telescopes. 

Among the large telescopes currently in various stages 
of dev elopment are the Giant Magellan Telescope (GMT, 
iJohns et all 2004. 20- 25m), C alifornia Extre mely Large Tele- 
scope (CELT, lNelsonll200 0b. 30m), Euro50 dAndersen et all 
[ 20041 50ml Over whelmingly Large Telescope (OWL, 
Dierickx et al. 2004, 100m), Japanese Extremely Large Tele- 
scope (live et alJl2004t 30m) ? and t he Chinese Future Giant 
Telescope (Ding-Oiang et al. 2004, 30m). The planned ca- 
pabilities of these telescopes are impressive. For instance, 
the GMT design calls for a resolution of #glao ~ 0.15" with 
ground layer adaptive optics (GLAO) and #fao ~ 0.007" with 
fully adaptive optics (FAO), corresponding to physical scales 
at redshift z = 3 of /?glao ~ 115 kpc and Rp\o ks 54 pc, re- 
spectively. The remnant systems in our simulations appro- 



priate for z = 3 progenitors have effective radii ranging from 
R e ff m 0.3- 10 physical kpc, making central stellar velocity 
dispersion measurements feasible in principle, and likely pos- 
sible even with GLAO observations for the most massive sys- 
tems. Infrared spectroscopy of rest-frame UV stellar lines 
could be used to find central stellar velocity dispersions of 
spheroids and, combined with the simulation results to es- 
timate Mbh via an inversion of the calculated Mbh-c rela- 
tion for a given redshift, the cosmic black hole po pulation 
can be constrained (see, e.g. lYu & Tremaindl2002l) . Future 
lensing observations with the Square Kilometer A rray to mea- 
sure black hole masses at intermediate redshifts ( Rusin et al. 
2005) may also be useful in this regard. In addition, if SMBH 
feedback is imp ortant to determining the properties of norma l 
galaxies (e.g. [Springel et al. 2005b; Robertson et al. 2005), 
such observations could prove to be crucial input for theo- 
ries of galaxy formation during the evolutionary epoch when 
the most massive black holes are formed. 

5. DISCUSSION 

Our simulations demonstrate that the merging of disk galax- 
ies whose properties are appropriate for redshifts z = - 6 
should produce an Mbh-c relation with the observed power- 
law scaling M B h oc cr 4 . The Mbh-c relation retains this scal- 
ing at redshifts z = - 6, even while the merger progenitors 
range widely in their structural properties with redshift. A 
possible weak evolution in the normalization of the M^-a 
relation is found, which may result from an increasing veloc- 
ity dispersion owing to the redshift evolution of the potential 
wells of galaxies. Simple arguments for how such evolution 
would affect the velocity dispersion provide a similar redshift 
scaling at redshifts z = 2 - 6 to those calculated from the sim- 
ulations. As these redshift scalings systematically place high- 
redshift systems below the z = M B h-ct relation, subsequent 
accretion processes may further reduce the scatter in the ob- 
served Mbh-c relation over time. 

Our results also provide guidance to semi-analytical mod- 
els about the behavior of the Mrh-ct relation. We suggest t hat 
semi-analytical mo dels in the vein of Volonte ri et al J (12003) or 
Menci *etal] ( 1200 3 ) could be altered to explicitly maintain the 
local M B H-Mbui g e relation at redshifts z < 6, produce the same 
Mbh-c scaling at redshifts z < 6, and reproduce the local 
Msn-a normalization. Simultaneously such models should 
tie morphological information, such as spheroid growth, to 
the growth of super massive bl ack holes (for a related discus- 
sion, see |Bromlev et al. 2004). While some semi-analytical 
models include prescriptions for black hole accretio n or feed- 
back (e.g. lMenci et al.ll2003t IWvithe & Loeb1l2003l) . it is dif- 
ficult to graft all the effects of feedback onto these models 
without a detailed treatment of the complex interplay between 
gas inflows, starbursts, and black hole activity. These mod- 
els could also be altered to more fully capture the black hole 
mass-dependent effects of feedback on galaxies and their en- 
vironments (see also, e.g. jHopkins et all 20051. 

Overall, our simulations provide more general conclusions 
about the cosmological evolution of structure. The spheroidal 
components of galaxies are likely formed in the same process 
that sets the M^w-cr relation and therefore the demograph- 
ics of black holes traces the demogra phics of spheroids ove r 
cosmic time, as is inferred locally (lYu & Tremaind 12002). 
With the possible shared origin of spheroids and SMBHs 
in mind, we note that observations suggest the stellar pop- 
ulatio ns of cluster ellipticals f orm at redshifts z > 2 (e.g. 
van Dokkum & Stanford! l2003t iGebhardtet alJl2003l) . This 
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era coincides with the epoch z ~ 3 during which the most 
massive SMBHs are formed , as inferred from the evolution in 
the qu asar number density ( Schmi dTet alJll995l iBovle etafl 
2000). Future comparisons with high-redshift AGN and mas- 
sive ellipticals could help provide our modeling of individual 
mergers with a broader cosmological context. By connecting 
with the evolution of the observe quasar luminosity function 
th rough the P(zA IorM rh) supplied by calculations in the vein 
of Hopkins et al. (2005e), our model may account for cosmo- 
logical considerations such as evolution in either the mass- 
dependent galaxy merger rate, the cold gas fraction, or the 
global star formation rate. 

While uncertainties remain in our method, an analysis sim- 
ilar to that presented here could be performed for simulations 
with different treatments of the relevant physics of black hole 
accretion and feedback, star formation, or ISM processes. As 
the observational data improve, the comparisons presented in 
this paper could be used to constrain the physical prescrip- 
tions used in simulations of galaxy formation. For instance, 
our knowledge of possible redshift-dependence in the ther- 
mal coupling rytheim of black hole feedback to the gas is lim- 
ited. The normalization of the Mbh-c relation depends on 
f/theim for a given accretion model, and a strong increase in 
Vtheim would introduce a correspondingly larger scatter in the 
Msft-a relation at z = 0. A decrease in the thermal coupling 
with redshift would decrease the scatter at z = 0, but in prin- 
ciple this effect could be counteracted if the dominant mode 
of black hole accretion is more efficient at higher redshifts. 
Better observational samples will help inform our modeling 
of these physics by improving constraints from the statistical 
moments of the z = M^-a relation. 

While our method for comparing evolution in the Mbh-c 
relation with observational constraints involves adopting a 
probability distribution function P(z/\ logMeH) relating black 
hole mass with a characteristic formation redshift, the gen- 
eral formalism does not depend on the particular choice 
of P(z/\ logMBH) or the corresponding interpretation of the 
quasar luminosity function. A similar analysis may be re- 
peated for other theoretical frameworks connecting galaxy 
formation to supermassive black hole growth. 

The largest final black holes masses we consider are Mbh ~ 
3 x 10 9 M Q , which are produced in the mergers of the very 
largest galaxies in our simulations. In principle, the rarest 
gas-rich mergers at high-redshifts could produce extraordi- 
narily massive black holes. The upper mass limit for black 
holes in the simulations presented here only reflects the fi- 
nal potential well depth of the largest galaxies we chose to 
merge as the M^-a relation generated by these mergers is 
a self-regulated process. The black holes grow in size to ac- 
commodate the growing potential fomented by star formation 
in the highest density regions of the galaxy during the height 
of the merger. As the energy input from the black hole into the 
surrounding gas saturates the ability of the gas contained by 
the potential to cool, a strong outflow is driven from the cen- 
tral regions and the growth of the black hole and spheroid are 
truncated. There is no indication in our simulations that this 
process cannot continue in suitably massive gas-rich galaxy 
mergers to produce ~ 10 I0 M Q black holes, and future work 
to determine the possibility of forming such black holes in 
galaxies mergers is planned. 

However, evide nce from the red sequence of the galaxy 
color bimod alitv JStrateva et alJ I200U iBlanton et alJj2003t 
iBaldrv et all 12004 iBell et al.1 12004 iFaber et alJ 120051) sug- 
gests that the most massive galaxies are produced through 



the dry merging o f spheroids at z < 1 (iBell et al.l 120051 
Ivan DokkunJl2005S) . In this picture, 1O 1O M black holes 
would likely have to grow in mass partially through binary 
black hole mergers. Whether black hole coalescence can 
occur in collisionless systems on a reasonable times cale is 
still an unresolved theoretical proble m (e.g., Begelma n et alJ 
ll980tlMilosavIievic & MerrittlhoOll) . Cosmological simula- 
tions that include prescriptions for black hole growth may 
naturally account for the formation of the most massive black 
holes and spheroids, but even this approach will likely expe- 
rience limitations owing to simulated volume size as the most 
massive galaxies are extremely rare. Semi-analytical mod- 
els that include prescriptions for black hole feedback provide 
another promising alternative, as such calculations can take 
advant age of N-body simulations like the Millennium Run 
(Springel et al. 2005cJ) that simulate cosmological volumes 
with 10 10 dark matter particles. However, the semi-analytic 
approach may ultimately need to be calibrated against sim- 
ulations like those presented here if there is any hope of re- 
alistically capturing the complicated hydrodynamical process 
that generates the Mbh-ct relation. The existence and possi- 
ble origin of the most massive SMBHs therefore promises to 
be an interesting avenue of future research. 

Lastly, the robust nature of the M^-a relation produced 
by our simulations suggests that the merging of galaxies over 
cosmic time will produce a relation with small intrinsic scat- 
ter, in accord with observations. The weak evolution of the 
calculated Mbh-c normalization implies that the small ob- 
served scatter in the Mr h-ct reflects the robust phys ical ori- 
gin of the relation (e.g. Haehnelt & Kauffmann 2000). While 
the dis persion Am bh = 0.25-0.3 measured bv lTremaine et all 
( 2002) is estimated to be intrinsic to the M^-a relation, if the 
observational errors have been significantly underestimated 
then our physical model may need improvement. We note that 
prelim inary evidence of a second parameter in the M^-a re- 
lation ( Siopis et al. 2004) may be connected to statistical vari- 
ance introduced by a redshift-dependent normalization, which 
may take the form of either a correlation with galaxy mor- 
phology or environment. Further comparisons between the 
observational data and simulations could help clarify these 
possibilities, especially in a cosmological setting. The red- 
shift evolution of the Faber-Jackson relation measured in the 
remnants may bear on observations of the funda mental plane 
JDressler etal.ll 19871 iDjorgovski & DavislH987t) and will be 
explored in forthcoming work. 

6. SUMMARY 

We calculate the M^-a relation produced by the merging 
of disk galaxies appropriate for redshifts z = 0,2,3, and 6 us- 
ing a large set of 1 12 hydrodynamic simulations that include 
the effects of feedback from accreting supermassive black 
holes and supernovae. We develop a method for comparing 
the Mbh-c relation produced by simulations with existing ob- 
servational constraints through the statistical mome nts of the 
measu red relation using the formalism presented bv lYu & Lul 
(2004), and find the simulations to be in agr eement w ith the 
dispersion Ai og M BH = 0.25-0.3 reported by Tremai ne et alJ 
(2002). The conclusions of our work are: 

• The simulations suggest that the normalization of the 
Mm\-(j relation may experience a weak evolution 
with redshift while the slope is consistent with be- 
ing constant near the locally observed value of ~ 4 
(Ferra rese & Merrittll2000l iGebhardt etafl feOOO). As 



12 



Robertson et al. 



our simulations cover a wide range of virial mass, gas 
fraction, ISM pressurization, dark matter halo concen- 
tration, and surface mass density at each redshift, our 
results suggest the physics that set the M^-a scal- 
ing are remarkably insensitive to the properties of pro- 
genitor galaxies in the scenario where the Mbh-c rela- 
tion is created through coupled black hole and spheroid 
growth in galaxy mergers. 

• We demonstrate how to compare theoretical models for 
the M BH -cr relation that predict redshift evolution in ei- 
ther the slope or normalization to the observed M^-a 
relation by calculating the z = dispersion Ai og M B H °f 
the relation given a functional form for the evolution 
and a probability distribution function P{zf\\ogM-Rn) 
that relates black hole mass to a characteristic forma- 
tion redsh ift. Adopting a model for P{zf\\ogM^n) 
from lHopkins et ail (2005e). we calculate the allowed 
strength of redshift evolution given power-law, red- 
shift evolution in either normalization or slope con- 
sistent with the observed scatter at z = 0. The mod- 
est normalization evolution calculated by our simula- 
tions produces a scatter similar to that observed. We 
demonstrate how redshift evolution in the M^-a re- 
lation introduces mass-dependent dispersion and skew- 
ness about the average relation, and show how these 
statistical moments can be used to differentiate theoret- 
ical models for the M^-a relation that predict distinct 
redshift evolutions. 

• We relate our simulations to observations of the M^-a 
relation in AGN at redshifts z > and, given the 
observational and theoretical uncertainties, find the 
properties of our inactive remnants to be more con- 
sistent with t he redshift-in dependent M^-a relation 
measured by Shields et alj (|2( | )( [ >3h t han th e claims of 
iTreu et all J2004I) and IWalter et al.1 J2004I) that super- 
massive black holes may form earlier than galactic 
spheroids. We show that if the off set <51ogg = -0.16 
in the M^-a relation measured bv ITreu et al J |2004) 
for AGN at z = 0.37 is attributed to a power-law evolu- 
tion in the normalization of the velocity dispersion as 
a(z) oc a(z = 0)(1 +zr, the resultant scatter expected in 



the z = Mbh-c relation would be much larger than is 
observed. Additional comparisons between the active 
phase of our merger simulations and the observations 
of active galaxies are planned in future work. 

• Our simulations provide additional motivation for fu- 
ture observations of high-redshift normal galaxies us- 
ing the next-generation of large telescopes for the pur- 
pose of constraining the evolution of the cosmic black 
hole population with time. The planned 20- 30m-class 
telescopes with infrared spectrographs should be able to 
measure central stellar velocity dispersions in galaxies 
to z ~ 3. Our simulations can be used to connect those 
measurements to estimates for the supermassive black 
hole mass function and the distribution of black hole 
formation redshifts P(zj | logMeH)- A complete census 
of the cosmic black hole population will require mea- 
surements of the spheroid population at high redshifts 
and will serve as important input into structure forma- 
tion theories that account for the effects of black hole 
feedback on normal galaxies. 
With a method for constraining theories of the M^-a rela- 
tion with the statistical moments of the local M-^-a relation 
based on their predictions for redshift evolution in the rela- 
tion, future opportunities for direct comparisons between ob- 
servations and theoretical models should be plentiful. As the 
statistics of the observational results improve, the theoretical 
models will be able to adjust their physical prescriptions to 
better accommodate empirical constraints. Research into the 
Mbh-c relation should then continue to prove an important 
area of interest for future efforts to construct a more complete 
theory of galaxy formation. 
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